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ABSTRACT 

Results  are  presented  of  sphere  drag  measurements  made  in  the  free-molecular  and 
transitional  flow  regimes.  The  drag  data  were  obtained  using  a  drag  balance  and  the 
fre e-flight  technique.  Conditions  for  which  measurements  were  made  are  within  the 
ranges  1  <  Kn.  <  32,  2.9  <  M.  <  11.2,  0.8  <  Tw/T.  <  16.0,  and  2.5  <  S»  <  3.2.  An 
analysis  of  the  drag  on  a  sphere  in  rarefied  flow  is  also  presented,  and  the  results  are 
compared  with  experimental  data  and  other  theories.  The  analytical  results  well  predict 
the  trend  of  the  experimental  data  in  the  transition  regime  and  remain  valid  at  Knudsen 
numbers  for  which  previous  theories  are  not  applicable. 
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SECTION  I 
INTRODUCTION 

Investigations  of  atmospheric  winds,  density,  pressure,  and  temperature  using 
spherical  satellites  and  falling  balloons  (Ref.  1)  have  stimulated  interest  in  sphere  drag 
from  subsonic  to  hypersonic  flow  and  from  the  continuum  to  the  free-molecular  regime. 
Besides  the  dependency  of  drag  on  the  Mach  number,  or  speed  ratio,  drag  in  the 
free-molecular  and  transitional  flow  regimes  depends  significantly  on  the  sphere  wall  and 
free-stream  gas  temperatures.  Also,  there  is  the  possibility  that  various  materials  will  have 
different  accommodation  coefficients  at  satellite  velocities  and  thereby  influence  the  drag. 
An  example  of  a  body  which  may  have  all  of  these  factors  influencing  the  drag 
simultaneously  would  be  a  satellite  designed  with  surface  materials  of  different 
emissivities  to  provide  thermal  control.  An  accurate  measurement  of  the  drag  of  such  a 
body  is  very  difficult  because  (1)  in  order  to  simulate  high  Knudsen  number  (Kn.  = 
\J d)  flow  in  present  low  density  test  facilities,  the  model  must  be  relatively  small  and 
this  makes  actual  simulation  of  the  model  difficult,  and  (2)  the  actual  surface  condition 
of  the  materials  is  difficult  to  reproduce  even  if  the  same  materials  and  paints  are  used, 
because  of  condensed  gases,  etc.,  on  the  surface  of  the  model  which  may  not  be  present 
in  space.  The  effect  of  wall  temperature  on  drag  can  be  measured  in  wind  tunnel 
experiments,  and  the  effect  of  surface  condition  on  drag  is  usually  not  as  large  and  can 
be  investigated  better  using  a  molecular  beam.  Therefore,  the  effect  of  surface  condition 
on  drag  will  not  be  discussed  other  than  the  presentation  of  data  taken  using  a  sphere 
with  a  black  and  then  a  gold  surface. 

The  purpose  of  this  report  is  to  present  the  results  of  sphere  drag  measurements 
made  in  the  free-molecular  and  transitional  flow  regimes.  A  theoretical  analysis  of  the 
drag  on  a  sphere  in  the  transition  regime  is  also  presented,  and  the  theoretical  and 
experimental  results  are  compared.  This  analysis  predicts  the  free-molecular  drag  value  as 
Kn_  -*bo  and  the  Newtonian  drag  value  as  Kn»  0.  The  present  results  are  compared  with 
other  analytical  and  experimental  investigations  (Refs.  2  through  9). 

SECTION  II 

EXPERIMENTAL  PROCEDURE 

The  experimental  data  were  taken  in  the  Aerospace  Research  Chamber  (ARC) 
(8V)  of  the  Aerospace  Environmental  Facility  (AEF)  at  AEDC.  This  chamber  is  10  ft  in 
diameter,  20  ft  in  length,  and  uses  cryogenic  surfaces  to  pump  the  gas  from  the  wind 
tunnel  nozzle.  The  principal  pumping  capacity  of  the  ARC  (8V)  chamberwas  provided  by 
620  ft2  of  77°  K  liquid-nitrogen  cryosurfaces  and  240  ft2  of  15  to  20°  K  gaseous-helium 
cryosurfaces.  A  total  of  4-kw  gaseous-helium  refrigeration  capacity  was  used.  This 
arrangement  permitted  the  continuous  pumping  of  nitrogen  or  argon  at  a  mass  flow  rate 
of  10  gm/sec  at  a  chamber  pressure  of  about  10'5  mm  Hg.  The  pumping  requirements  for 
this  test  were  less  than  10  gm/sec. 
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Drag  measurements  were  made  using  a  balance  and  a  tree-flight  technique.  The 
design,  calibration,  and  operation  of  the  equipment  and  instrumentation  necessary  for 
these  two  methods  of  measuring  drag  are  discussed  in  detail  in  Ref.  10.  The  data 
reduction  techniques  are  discussed  here. 

2.1  WIND  TUNNEL  NOZZLES 

Two  low  density  nozzles,  designated  as  the  Mach  number  3  (M3)  and  Mach 
number  6  (M6)  nozzles,  were  used  to  provide  the  flow  conditions.  The  wall  cooling  on 
the  M6  nozzle  was  recently  modified  by  replacing  the  aluminum  tubing  through  which 
liquid  nitrogen  was  forced,  by  a  complete  liquid-nitrogen  jacket  similar  to  the  M3  nozzle. 
The  result  was  to  produce  a  constant  89°K  wall  temperature  which  was  a  100°K  decrease 
in  wall  temperature  at  the  nozzle  exit.  This  reduction  in  wall  temperature  reduced  the 
boundary-layer  growth  on  the  nozzle  wall  and  also  improved  the  repeatability  of  the  flow 
in  the  nozzle  since  the  wall  temperature  remained  constant  with  time.  Some  results  on 
the  flow  conditions  in  the  M3  nozzle  using  nitrogen  or  argon  as  the  test  gas  are  given  in 
Ref.  11.  Boundary-layer  profiles  and  axial-flow  conditions  for  the  M3  nozzle,  and  more 
results  on  the  flow  conditions  in  the  M3  and  M6  nozzles,  for  nitrogen  only,  are  given  in 
Refs.  12  and  13.  A  summary  of  the  operating  ranges  of  these  nozzles  is  given  in  Fig.  1 
(Appendix  I). 

2.2  NOZZLE  CALIBRATION 

The  total  pressure,  P0,  was  obtained  using  an  MKS  Baratron  ®  (variable 
capacitance  transducer),  which  measured  the  average  of  the  pressure  at  two  locations  in 
the  reservoir  as  shown  in  Fig.  2.  The  total  temperature,  T0 ,  was  measured  using  a  total 
temperature  probe  located  on  the  centerline  of  the  reservoir  (Fig.  2).  The  flow  was 
calibrated  using  a  l-in.-diam,  10-deg  internally  chamfered  pitot  probe.  A  Baratron  was 
used  to  measure  the  pitot  pressure.  The  pitot  probe  was  mounted  on  a  scanner 
mechanism  as  shown  in  Fig.  2.  The  movement  of  the  scanner,  which  was  in  the  plane  of 
the  figure,  permitted  the  measurement  of  pitot  pressures  from  5  in.  inside  to  20  in. 
downstream  of  the  nozzle  exit  planes  and  from  S  in.  below  to  20  in.  above  the  nozzle 
centerlines.  The  vertical  movement  was  sufficient  to  traverse  the  exit  radii  of  the  two 
nozzles. 


For  the  rarefied  flow  conditions  of  interest,  the  viscous  corrections  to  the  pitot 
pressure  measurements  were  significant.  The  magnitude  of  this  effect  on  Mach  number 
and  dynamic  pressure,  q^,,  is  illustrated  in  Fig.  3  for  the  M6  nozzle.  Since  the  drag 
coefficient  is  inversely  proportional  to  q_,  it  is  obvious  from  Fig.  3  that  serious  error  is 
introduced  if  the  viscous  correction  is  not  known  accurately.  The  viscous  correction  data 
used  to  correct  the  pitot  pressure  measurements  were  obtained  in  a  previous  test 
program,  and  these  data  are  given  in  Fig.  4  with  data  from  other  sources  (Refs.  14 
through  17). 

A  problem  associated  with  low  density  wind  tunnel  testing  is  that  of 
boundary-layer  merging,  i.e.,  the  boundary  layer  completely  filling  the  nozzle.  It  has  been 
previously  assumed  that  merging  in  the  M6  nozzle  occurred  at  the  point  of  minimum 
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Mach  number  (Fig.  3).This  assumption  has  since  been  found  to  be  invalid  for  the 
following  reason.  It  was  observed  that  an  abrupt  increase  in  the  drag  coefficient  occurred 
at  values  of  Kru  where  the  drag  coefficient  should  be  free  molecular  or  at  least  be 
asymptotically  approaching  the  free-molecular  limit.  In  order  to  explain  the  unexpected 
increase  in  Cp,  an  investigation  was  made  to  determine  whether  or  not  the  flow  was 
merged,  because  if  it  was,  the  pitot  pressure  measurements  would  be  incorrect  and  the 
free-stream  properties,  q_,  T„,  etc.,  calculated  using  P0,  T0)  and  P0',  would  also  be 
incorrect.  The  problem  of  merged  flow  was  investigated  using  the  fully  viscous 
slender-channel  program  developed  by  Rae  (Ref.  18).  It  was  found  from  the  numerical 
solutions  that  a  significant  decrease  in  the  local  total  enthalpy,  H,  occurred  at  the  M6 
nozzle  exit  centerline  at  about  the  same  point,  i.e.,  value  of  P0  for  T0  =  constant,  at 
which  the  increase  in  Cp  occurred  (Fig.  5).  The  decrease  in  Hq,  e  indicated  the  flow  was 
not  isentropic,  and  hence  the  q„  used  to  determine  Cp  for  P0  <  1  mm  Hg  was 
incorrect.  The  value  of  nozzle  reservoir  Reynolds  number,  Reo  r*,  at  which  merging 
occurred  for  this  value  of  the  nozzle  wall  to  reservoir  temperature  ratio  (which  was  0.1 1), 
was  Re0>r*  =  900.  The  throat  radius  for  this  nozzle  was  r*  =  1.54  in.  It  is  interesting  to 
note  that  merging  actually  occurred  for  the  conditions  in  Fig.  3  at  PQ  «  1  mm  Hg,  which 
was  a  significantly  larger  value  of  P0  than  the  value  at  which  the  minimum  Mach  number 
occurred.  Numerical  solutions  were  also  obtained  for  the  flow  in  the  M3  nozzle;  however, 
data  are  presented  only  for  isentropic  flow  conditions. 

2.3  DATA  REDUCTION 

2.3.1  Drag  Balance 

Data  reduction  of  the  forces  measured  using  the  drag  balance  was  straightforward. 
After  force  measurements  were  made  on  the  sphere  mounted  on  the  sting  at  various  flow 
conditions,  the  sphere  was  removed,  and  a  second  sphere,  independently  supported,  was 
positioned  slightly  upstream  of  the  sting  (Fig.  6).  The  tare  forces  on  the  sting  were 
measured  using  this  arrangement  for  the  same  range  of  flow  conditions  for  which  total 
force  measurements  were  made.  An  example  of  the  tare  force  and  total  drag  force 
measurements  is  given  in  Fig.  6.  For  this  0.25-in. -diam  sphere,  the  tare  force  was  from  15 
to  20  percent  of  the  total  drag  force.  If  the  sphere  size  was  reduced,  the  percent  of  total 
force  which  was  attributable  to  tare  increased.  The  difference  between  tare  and  total 
force  measurements  is  the  sphere  drag,  D.  The  drag  coefficient,  Cp,  was  obtained  from 
the  definition 


CD 


D 

!/2  PJ £  A 


D 

‘UoA 


(1) 


A  second  sphere,  identical  to  the  one  used  on  the  balance  sting  to  measure  the 
drag,  was  used  to  determine  the  sphere  wall  temperature.  This  second  sphere  was 
supported  by  the  leads  to  a  thermocouple  which  was  used  to  measure  the  sphere 
temperature.  It  was  attached  to  the  same  support  as  the  balance  and  was  simultaneously 
subjected  to  the  same  flow  conditions  as  the  sphere  mounted  on  the  balance  sting.  It  was 
assumed  that  the  temperature  of  the  sphere  was  uniform,  and  therefore  the  measured 
temperature  was  taken  as  the  sphere  wall  temperature. 
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2.3.2  Free  Flight 

Free-flight  drag  data  were  obtained  from  photographs  of  the  trajectories  of 
spheres  (Fig.  7)  which  were  dropped  through  the  nozzle  test  section.  The  photographs 
were  multiple  exposures  of  spheres  illuminated  at  constant  time  intervals  by  a  strobe  light 
operating  at  the  rate  of  120  flashes/sec.  From  one  to  three  aluminum  and/or  magnesium 
spheres  were  dropped  simultaneously.  The  spheres  were  initially  located  in  the  drop 
mechanism  a  minimum  of  ten  diameters  apart  in  a  plane  perpendicular  to  the  plane  of 
the  photograph  in  Fig.  7.  Immediately  after  the  spheres  were  dropped,  the  minimum 
separation  distance  increased,  and  it  was  assumed  that  no  interaction  among  the  spheres 
occurred.  Sphere  diameters  of  1/32,  1/16,  and  1/8  in.  were  used.  The  sphere  wall 
temperatures  were  determined  by  measuring  the  temperature  of  the  drop  mechanism.  The 
change  in  sphere  temperature  for  the  period  of  time  required  for  the  spheres  to  traverse 
the  test  section  was  negligible. 

The  drag  force  on  a  sphere  was  determined  by  applying  Newton’s  second  law  of 
motion  to  the  sphere  in  the  direction  of  the  free-stream  velocity.  The  distance,  x,  which 
a  constant  mass  body  would  move  in  a  constant  force  field  is  a  quadratic  equation  in 
time,  t,  i.e., 

X  =  cit2  +  c2t  +  c3 

Using  Newton’s  law 

and  substituting  Eq.  (2)  into  Eq.  (3)  and  differentiating  gives 

D  =  2mci  (4) 

where  m  is  the  mass  of  the  sphere.  Therefore,  if  ci  can  be  determined,  the  drag  can  be 
calculated.  The  constant  ci  was  determined  by  fitting  a  second-order  polynomial  through 
the  data  of  a  plot  of  x  versus  t  (Fig.  8).  The  value  of  x  for  each  t  was  determined  by 
measuring  the  axial  displacement  of  the  sphere  from  point  to  point  in  Fig.  7  at  constant 
(known)  intervals  of  time.  These  measurements  were  made  from  the  photographic 
negatives  using  a  film  reader. 


(2) 

(3) 


SECTION  III 

TRANSITION  DRAG  ANALYSIS 

A  simple  analysis  is  presented  for  the  drag  of  spheres  in  the  transitional  flow 
regime.  The  mathematical  model  is  based  on  the  following  four  assumptions:  (1)  The 
flow  Mach  number  is  sufficiently  high  that  the  number  flux  of  molecules  which  has  a 
chance  of  colliding  with  the  sphere  is  given  by  N„  =  n»U„,  (2)  all  molecules  are  emitted 
normal  to  the  sphere  surface,  (3)  each  emitted  molecule  experiences  only  one  collision 
with  an  oncoming  free-stream  molecule,  after  which  neither  molecule  is  considered,  and 
(4)  all  collisions  occur  at  a  distance  Xw  from  the  sphere  surface  (Fig.  9). 
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There  are  certain  restrictions  and  conditions  which  need  to  be  considered  because 
of  the  basic  assumptions.  The  first  assumption  requires  that  the  random  molecular 
velocity  be  less  than  the  mean  flow  velocity.  Physically  this  means  that  the  molecules 
tend  to  move  in  straight  lines  parallel  to  the  mean  flow  direction,  and  therefore,  the  net 
number  which  cross  the  boundary  of  the  imaginary  cylinder  swept  out  by  the  sphere  per 
unit  time,  i.e.,  irr2Um,  is  negligible.  This  is  not  a  particularly  stringent  assumption,  as 
indicated  by  the  results  given  in  Fig.  10  for  the  incident  flux  in  free-molecular  flow  (see 
Ref.  19  for  a  derivation  of  the  incident  flux).  From  Fig.  10  the  approximation  Nj  *  N„ 
=  n_U„  is  good  for  S*  >  1  or  (2RTW)^  <  U„. 

The  second  assumption  permits  one  to  neglect  the  distribution  function  of  the 
molecules  emitted  from  the  sphere  surface,  and  it  is  justified  only  by  the  final  results. 

The  third  and  fourth  assumptions,  together,  imply  that  the  sphere  is  in 
free-molecular  flow  with  respect  to  molecules  which  survive  the  collision  surface  (Fig.  9) 
since  there  are  no  collisions  between  free-stream  and  emitted  molecules  in  the  region 
between  the  collision  surface  and  the  sphere.  Therefore,  by  determining  the  number  flux 
of  surviving  molecules,  nsU„,  the  drag  on  a  sphere  can  be  determined  from 

D  Cd  ns 

Dfm  CDfm  noo  ^ 

since  the  drag  on  the  sphere  is  the  same  as  the  free-molecular  drag  with  a  free-stream 
number  density  of  ns.  That  is,  for  the  same  free-stream  and  sphere  wall  conditions,  the 
drag  is  proportional  to  the  free-stream  number  density,  or,  for  the  case  of  the  surviving 
molecules,  to  the  effective  free-stream  number  density,  n,. 

Before  solving  for  nsU„,  it  is  interesting  to  note  an  error  introduced  by 
considering  the  ratio  CD  /CDfm  in  Eq.  (5).  It  was  assumed  in  writing  Eq.  (5)  that  the 
free-stream  conditions  of  the  surviving  molecules,  other  than  number  density,  did  not 
change,  and  also  that  the  sphere  wall  conditions  did  not  change.  The  assumption  of 
constant  free-stream  conditions  is  valid  because  of  the  four  basic  assumptions,  but  the 
one  concerning  the  wall  conditions  is  not.  This  is  because  fewer  molecules  are  available  to 
transfer  heat  between  the  gas  and  the  sphere  because  of  the  “shielding”  effect  produced 
by  the  collisions,  and  therefore,  Tw  is  not  the  same  for  the  shielded  and  unshielded 
sphere.  However,  for  a  diffusely  reflecting  wall,  which  was  the  model  assumed  in 
calculating  Cofm  for  the  experimental  data,  the  effect  of  Tw  decreases  as  S«,  increases 
(Fig.  11),  and  it  is  assumed  that  the  effect  of  the  difference  between  wall  temperatures 
on  Cq  is  negligible.  Of  course  for  a  specularly  reflecting  wall,  Cp  is  independent  of  Tw, 
and  no  error  is  introduced. 


In  the  solution  for  nsU„  the  mean  distance  between  collisions  of  emitted 
molecules  with  oncoming  molecules,  Aw ,  will  be  taken  as  a  function  of  ip  as  given  by 
(Refs.  2  and  20).. 


Kf  = 


Vw 

noo  [(Uoc  +  vw  cos  ^)2  +  (vw  sin  ^)2]^ 


(6) 
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In  terms  of  the  mean  free  path  of  the  molecules  in 
Eq.  (6)  can  be  written 

Aw  yj2  vw 

Aeo  [Uoo  *■  2U00VW  COS  t lj 


the  free  stream.  A,  =  Xlisflita2^), 


+  V, 


,2]* 


(7) 


Since  the  number  flux  which  has  a  chance  of  colliding  with  the  sphere  is  n„U„, 
then  only  those  collisions  which  occur  at  the  “effective  collision  surface”  (which  is  that 
part  of  the  collision  surface  within  the  cylindrical  volume  swept  out  by  the  sphere,  Fig. 
9)  will  affect  the  surviving  number  _flux  which  reaches  the  sphere.  The  number  of 
molecules  per  unit  time  which  survive,  Ns,  can  be  expressed  as 

Ns  =  No*  -Nc  (8) 

where 

Nm  =  n^Uoenr2  (9) 


is  the  number  of  free-stream  molecules  per  unit  time  which  would  collide  with  the  sphere 
if  no  collisions  between  emitted  and  free-stream  molecules  occurred,  and  Nc  is  the 
number  of  molecules  per  unit  time  which  are  emitted  toward  the  effective  collision 
surface  and  hence  suffer  collisions  with  free-stream  molecules.  Therefore,  Nc  is  the 
number  per  unit  time  which  are  removed  from  N„.  Since  the  molecules  were  assumed  to 
be  emitted  normal  to  the  sphere  surface,  the  number  of  molecules  emitted  from  the 
sphere  surface  within  the  angle  0  is  equal  to  Nc  as  expressed  by 

0 

Nc  =  ngUooAp  =  nsUoo  /  2»ttz  sin  iji  cos  if/  di£r  (10) 

O 

where  ns  is  the  number  density  behind  the  collision  surface,  and  the  integral,  Ap ,  is  the 
projected  surface  area  of  the  sphere  within  0  (Fig.  9).  Integrating  Eq.  (10)  gives 


Expressing  Ns  as 


Nc  =  ngUooTJT 2  sin2  6 
Ns  =  nsUoe?7r2 


(in 

(12) 


which  is  the  total  number  per  unit  time  which  survive  to  collide  with  the  sphere,  and 
then  substituting  Eqs.  (9),  (11),  and  (12)  into  Eq.  (8)  gives 


n§  Uoo 


nooUoo  ns  Uo 


Solving  for  ns/n*one  has 

ns  l 

Uoo  1  +  sin2  Q 


(13) 

(14) 


Substituting  Eq.  (14)  into  Eq.  (5)  and  using  Ap/Afor  sin2  0  (Eqs.  (10)  and  (11))  gives 


cD _ i 

^D{m  1  +  sin2  Q 


(15) 
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The  angle  6  depends  on  Xw  since  6  is  the  angle  at  which  the  distance  between  the 
collision  surface  and  sphere  surface  in  the  Q  direction  is  equal  to  Xw  (Fig.  9).  From  the 
geometry  in  Fig.  9  one  can  write 

sin  6  -  - - -  (16) 

r  +  \iJj=Q 

Substituting  Eq.  (16)  into  Eq.  (15)  gives 


C  p 

CDfm 


(17) 


Using  Eq.  (7)  for  Xw .  the  term  X*  /r  becomes 


Aw  (Knoot  ft 

r 


2y2  Kn©c 


(18) 


where  Kn„  =  XJd.  Sw  =  UJ(2RTW  )'A ,  and  vw  =  (9ffRTw/8)vi  (Ref.  5).  Substituting  Eq. 
(18)  into  Eq.  (17)  gives  for  the  drag  coefficient 


C.p 

rDfm 


An  expression  for  0  used  in  Eq.  (19)  can  be  obtained  by  substituting  Eq.  (18) 
into  Eq.  (16)  to  get 


sin  0  - 


1 


2y  2  K  n  oo 

’( 4S'w  V 

i  *>  fiSw]cn-  0  ' 

1  s 

3\h  ) 

\  3  \JttI 

(20) 


Equation  (20)  is  a  transcendental  equation  which  must  be  solved  for  6  for  each  value  of 
Kn„,  and  Sw.  The  solution  for  6  which  satisfies  0  <  6  <  ir/2  (Fig.  9)  was  used  in  Eq.  (19) 
to  calculate  Cd/Cd,,,,  for  each  value  of  Kn„  and  Sw. 

SECTION  IV 

RESULTS  AND  DISCUSSION 


4.1  ANALYTICAL 

It  is  of  interest  to  examine  the  limits  of  Eq.  (19).  In  the  limit  Kn. -soothe  result 
for  the  drag  coefficient  is  Cd/Cd^  =  1.  which  is,  of  course,  expected,  In  the  limit 
Kn„-+  0  the  result  for  the  drag  coefficient  is  CD/CDfm  =  0.5.  This  result  can  be  reasoned 
physically  by  recalling  that  one  of  the  assumptions  was  that  only  one  collision  between 
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emitted  and  free-stream  molecules  was  considered.  Therefore,  when  the  mean  free  path 
X_  =  0  (and  hence  Xw  =  0,  see  Eq.  (18)),  only  one-half  the  molecules  which  have  a 
chance  of  colliding  with  the  sphere  could  actually  collide,  since  each  molecule  at  the 
sphere  surface  prevented  another  molecule  from  colliding  with  the  surface.  This  result  is 
the  same  as  the  result  obtained  from  the  Newtonian  theory  if  CDfm  is  based  on  specular 
reflection.  Newtonian  theory  assumes  that  all  the  molecules  directly  upstream  of  the 
sphere  hit  the  sphere  but  that  they  impart  momentum  because  of  the  collision  only, 
since  they  leave  parallel  to  the  surface.  Whereas,  in  the  present  analysis,  only  one-half  the 
molecules  directly  upstream  of  the  sphere  hit  the  sphere,  but  they  impart  momentum 
from  both  the  collision  and  the  reflection.  This  is  the  reason  for  a  sensible  result  in  the 
continuum  limit. 

Some  results  from  Eq.  (19)  are  given  in  Fig.  12  where  0  was  obtained  from  the 
solution  to  Eq.  (20)  for  each  value  of  Kn.  and  Sw.  Notice  in  Fig.  12  that  Cd/Cbj,,, 
decreases  as  Sw  increases.  That  is,  for  constant  free-stream  conditions,  the  collision 
surface  approaches  the  sphere  surface  as  Sw  increases,  because  of  decreasing  Tw  and 
hence  decreasing  Xw ,  and  consequently  more  shielding  occurs,  and  the  drag  is  reduced. 

A  simplification  in  the  calculation  of  the  drag  coefficient  is  obtained  if  Xw  in  Eq. 

(17)  is  taken  as  a  function  of  Kn„  and  Sw  only,  i. e.,  if  6  in  Eq.  (18)  is  taken  as 
constant.  The  effect  of  the  6  term  in  Eq.  (18)  is  shown  in  Fig.  13  to  have  relatively  little 
effect  on  the  value  of  Xw,  or  Knw,  for  constant  Sw.  The  assumption  of  constant  9  in 
Eq.  (18)  eliminates  the  necessity  of  solving  Eq.  (20).  The  resulting  drag  coefficients  as 
calculated  using  constant  and  variable  Xw  are  illustrated  in  Fig.  14  for  Sw  =  1.  For  small 
values  of  Kn„  the  result  with  Xw(0  =  7t/2)  gives  slightly  better  agreement  with  variable 
Xw ,  whereas  for  large  values  of  Kn„  the  result  with  Xw  ( 9  =  0)  gives  better  agreement.  This 
result  is  expected  since  for  small  Kru,  9  *  r/2,  and  for  large  Kn..,  0*0  (Fig.  9).  In 
general,  for  all  Kn»,  the  result  with  Xw  (0  =  0)  gives  better  agreement  with  the  result  for 
variable  Xw  (Fig.  14).  It  should  be  pointed  out  that  the  assumption  of  constant  0  in  Eq. 

(18)  is  justified  only  because  Xw  is  a  weak  function  of  0  and  a  strong  function  of  Kn. 
and  Sw. 

4.2  EXPERIMENTAL 
4.2.1  Balance  Data 

Drag  measurements  were  made  using  the  balance  in  the  M3  and  M6  nozzle  flows. 
The  reservoir  (total)  temperature  was  800°K,  and  the  total  pressure  was  varied  between 
0.1  and  5.0  mm  Hg,  depending  on  the  nozzle  used,  to  produce  free-stream  Knudsen 
numbers  of  1  <  Kn«  <  5,  based  on  a  0.25-in.-diam  sphere. 

The  balance  proved  capable  of  measuring  the  change  in  sphere  drag  which  was 
caused  only  by  a  changing  wall  temperature.  The  wall  temperature  of  the  0.25-in.-diam 
sphere  used  to  take  the  data  in  Fig.  15  varied  between  253  and  314°K,  and  the 
free-stream  gas  static  temperature,  T„,  varied  between  74  and  62°K.  Taking  into  account 
the  appropriate  S„  and  Tw/T„  ratio  in  the  calculation  of  Cq  fm  ,  the  drag  data  normalized 
by  Cofm  have  less  deviation  than  the  Cq  data  as  shown  in  Fig.  15.  The  assumption  used 
for  the  calculation  of  Cd  f m  was  that  the  molecules  were  completely  accommodated  and 
reflected  diffusely  from  the  sphere  surface  (Fig.  1 1). 
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The  balance  drag  data  are  shown  in  Fig.  16  by  the  closed  symbols.  The  data 
taken  in  the  M6  nozzle,  i.e.,  7.0  <  <  7.7,  have  less  than  2-percent  scatter,  and  the 

data  taken  in  the  M3  nozzle  have  slightly  more  scatter.  The  difference,  if  any,  in  the  drag 
between  the  two  spheres  with  painted  surfaces,  one  coated  with  a  flatblack  lacquer  paint 
of  high  emissivity,  which  was  used  for  thermal  control  on  satellites,  and  the  other  coated 
with  gold,  could  not  be  detected. 

The  drag  coefficients  and  flow  conditions  for  the  balance  data  are  tabulated  in 
Table  I  (Appendix  II). 

4.2.2  Free-Flight  Data 

The  free-flight  drag  data  using  nitrogen  and  argon  as  the  test  gas  are  also  given  in 
Fig.  1 6.  The  scatter  in  the  free-flight  data  is  obviously  much  worse  than  the  scatter  in  the 
balance  data.  There  are  several  possible  explanations  for  the  scatter  in  the  free-flight  data. 
Among  them  are:  (1)  The  spheres  were  dropped  through  the  test  section  with  initial 
direction  normal  to  the  free  stream  rather  than  directly  upstream  as  is  usually  done  in 
free-flight  measurements,  and  therefore,  small  axial  displacements  of  the  spheres  were 
obtained  because  of  the  low  free-stream  dynamic  pressure  and  the  requirement  that  the 
spheres  be  at  least  1/32  in.  in  diameter  in  order  to  be  photographed,  (2)  the  strobe  light 
used  was  limited  to  120  flashes/sec,  and  this  restricted  the  number  of  exposures  of  the 
spheres  which  could  be  made  during  the  free  fall,  (3)  error  in  reading  the  sphere 
trajectories  from  the  photograph  negatives,  (4)  large  distance,  9  ft,  from  the  camera  and 
strobe  light  to  the  spheres,  and  (5)  the  portion  of  the  test  section  through  which  the 
spheres  fell  had  an  axial  Mach  number  gradient  of  about  0.05/in.,  and  the  flow 
conditions  used  to  reduce  the  sphere  data  corresponded  to  that  point  in  the  flow  where 
the  spheres  crossed  the  nozzle  axis  (the  maximum  axial  displacement  of  the  spheres  was 
about  6  in.).  Suggestions  to  eliminate  or  reduce  the  magnitude  of  these  effects  are  given 
in  Ref.  10. 

The  drag  coefficients  and  flow  conditions  for  the  free-flight  data  are  tabulated  in 
Table  I. 

4.3  ANALYTICAL  AND  EXPERIMENTAL  COMPARISONS 

Compared  with  the  balance  and  free-flight  drag  data  in  Fig.  16  is  the  analytical 
result  given  by  Eq.  (19)  for  Sw  =  2.7,  which  is  a  mean  value  of  the  wall  speed  ratio  for 
these  data.  From  this  comparison  the  analysis  appears  to  provide  a  reasonable  prediction 
of  the  drag  in  the  near  free-molecular  regime.  Certainly  the  agreement  between  theory 
and  experiment  is  good,  considering  the  simple  mathematical  model  used  to  analyze  the 
drag  in  the  transition  regime. 

In  order  to  determine  the  qualitative  trend  of  the  analysis  throughout  the 
transitional  flow  regime,  Eq.  (19)  is  compared  with  drag  data  from  various  sources  in  Fig. 
17  for  two  values  of  Sw  which  essentially  bound  the  experimental  data.  As  can  be  seen 
in  Fig.  17,  the  trend  of  the  experimental  data  in  the  transition  regime  is  relatively  well 
predicted.  However,  the  experimental  data  of  Fig.  1 7  do  not  always  vary  with  Sy,  as  the 
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analytical  results  predict.  Therefore,  the  solution  for  Sw  =  l'was  taken  to  compare  with 
some  data  in  Fig.  18  which  have  little  scatter.  It  was  found  that  the  solution  for  Sw  =  1 
provided  the  best  agreement  with  the  bulk  of  the  experimental  data.  Furthermore,  if  Xw 
is  taken  as  constant  with  respect  to  8,  then  the  curve  in  Fig.  18  can  be  rather  accurately 
approximated  by 


Cd 


Cd 


fm 


1 


1  1  2 
+  1.615  Knoo 


(21) 


which,  in  terms  of  M_,  Re„,  and  7>  can  be  written  as 


CD  i 


1  + 

1 

2 

,  .  2.02  Vy 

^  R  6  eo  ^ 

(22) 


A  comparison  is  made  in  Fig.  19  of  the  present  result  with  the  theories  of  Baker 
and  Charwat  (Ref.  2),  Rose  (Ref.  3),  and  Willis  (Ref.  4)  and  some  experimental  data  for 
near  1.6.  The  theories  of  Rose  and  Willis  depend  on  S„  and  are  calculated  only  for 
the  conditions  of  Sims’  data  (Ref.  7).  They  should  not  be  expected  to  compare  with  the 
other  data,  although  the  conditions  of  the  data  of  Ref.  6  are  close  to  those  of  Sims’  data 
(the  different  test  gases  used  should  be  noted,  however).  The  present  result  and  the 
theory  of  Baker  and  Charwat  do  not  depend  on  S„  and  therefore  are  applicable  to  all  the 
data  in  Fig.  19.  As  can  be  seen  in  Fig.  19,  the  theory'  of  Willis  and  the  present  result  well 
predict  these  drag  data  in  the  near  free-molecular  regime.  Also,  the  present  result 
continues  to  predict  these  drag  data  at  Knudsen  numbers  below  which  the  existing 
theories  shown  in  Fig.  19  are  not  applicable. 

SECTION  V 
CONCLUSIONS 

Sphere  drag  measurements  have  been  made  using  the  balance  and  free-flight 
techniques  within  the  ranges  1  <  Kic  <  32,  2.9  <  <  11.2,  0.8  <  Tw/T„  <  16.0,  and 

2.5  <  Sw  <  3. 2. The  balance  data  had  less  than  ±1 -percent  scatter,  but  the  free-flight  data 
had  as  much  as  ±  10-percent  scatter.  Drag  measurements  were  made  on  two  spheres  which 
had  different  surface  coatings,  one  painted  black  and  the  other  gold,  and  any  difference 
in  the  drag  coefficients  of  the  two  spheres  was  less  than  the  experimental  error  of  the 
balance  data. 

It  was  found  that  the  balance  used  was  sensitive  enough  to  measure  the  increase 
in  drag  on  a  0.25-in.-diam  sphere  which  was  caused  only  by  an  increasing  sphere  wall 
temperature.  Proper  accounting  of  Tw  in  calculating  Cofm  for  normalizing  CD  smoothed 
the  drag  data. 

The  sphere  drag  analysis  provided  a  reasonable  prediction  of  the  drag  coefficient 
throughout  the  transitional  flow  regime.  It  was  in  good  agreement  with  the  theory  of 
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Willis,  which  appeared  to  be  the  most  accurate  of  the  three  theories  with  which  it  was 
compared  in  the  near  free-molecular  regime  for  the  conditions  considered.  Furthermore, 
the  present  analysis  was  in  good  agreement  with  the  experimental  data  for  Knudsen 
numbers  below  which  the  previous  theories  were  not  valid. 
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Fig.  2  Schematic  of  Low  Density  Tunnel,  ARC  (8V) 
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Fig.  4  Pitot  Probe  Viscous  Correction  Data 
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Fig.  5  Increase  in  CD/CDfm  and  Decrease  in  Hq  /H_  for  the  M6 
Nozzle  (T0  =  800°  K) 
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Fig.  6  Tare  Force  and  Total  Drag  Measurements  Using  0.25-in.-diam  Sphere 
on  Drag  Balance  in  M6  Nozzle  Flow 
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Fig.  7  Falling  Spheres  Using  the  Free-Flight  Technique 
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Time,  Number  of  Strobes 


Fig.  8  Curve  Fit  of  the  Axial  Displacement  of  a  Sphere  versus  Time 
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Fig.  9  Assumed  Model  of  the  Rarefied  Flow  about  a  Sphere 
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Fig.  10  Incident  Flux  of  Molecules  on  a  Sphere  in  Free-Molecular  Flow 
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Fig.  13  Effect  of  0  on  Knw  /Kn„ 
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Fig.  15  Effect  of  Normalizing  CD  by  CDfm 
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Fig.  18  Analytical  Result  with  S*  °  1  and  Some  Experimental  Data  with  Little  Scatter 


AEDC-TR -70-32 


Eqs.  (19)  and  (20)  with  S  . -  1.6 
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TABLE  I 

SPHERE  DRAG  COEFFICIENTS  AND  FLOW  CONDITIONS 
a.  Balance  Data 


M. 

Re. 

Kn_ 

Tw/T. 

Sw 

Cp 

CD/CDfm 

Test  Gas 

6.98 

2.12 

4.92 

3.40 

3.16 

2.460 

1.012 

Nitrogen 

7.02 

2.30 

4.55 

3.49 

3.14 

2.437 

1.001 

7.06 

2.47 

4.26 

3.58 

3.12 

2.462 

1.011 

7.10 

2.68 

■3.95 

3.68 

3.09 

2.445 

1.003 

7.13 

2.82 

f  3.77 

3.77 

3.07 

2.450 

1.003 

7.19 

3.10 

•'  3.45 

3.89 

3.05 

2.460 

1.006 

7.25 

3.37 

3.21 

4.01 

3.03 

2.465 

1.008 

7.31 

3.59 

3.03 

4.15 

3.00 

2.455 

1.003 

7.41 

4.01 

2.75 

4.34 

2.97 

2.425 

0.990 

7.49 

4.75 

2.35 

4.51 

2.95 

2.380 

0.972 

7.57 

5.72 

1.97 

4.67 

2.93 

2.310 

0.942 

7.65 

7.40 

1.54 

•  4.87 

2.90 

2.210 

0.901 

7.70 

9.25 

1.24 

*  5.05 

2.86 

2.152 

0.875 

7.22 

3.20 

3.36 

4.08 

2.99 

2.480 

1.012 

7.40 

3-95 

2.79 

4.43 

2.94 

2.437 

0.993 

7.50 

4.81 

2.32 

4.70 

2.89 

2.418 

0.984 

7.57 

5.70 

1.98 

4.86 

2.87 

2.315 

0.942 

7.61 

6.51 

1.74 

5.00 

2.85 

2.265 

0.919 

7.65 

7.36 

1.55 

5.13 

2.82 

2.220 

0.900 

7.68 

8.30 

1.38 

5.25 

2.80 

2.172 

0.880 

7.70 

9.03 

1.27 

5.37 

2.78 

2.156 

0.871 

2.90 

1.19 

3.65 

0.80 

2.71 

2.615 

0.951 

2.92 

1.23 

3.55 

0.85 

2.65 

2.630 

0.955 

3.07 

2.16 

2.12 

0.90 

2.71 

2.405 

0.890 

3.19 

2.64 

1.80 

1.00 

2.67 

2.420 

0.890 

3.28 

3.23 

1.51 

1.10 

2.61 

2.250 

0.820 

3.34 

3.89 

1.28 

1.20 

2.55 

2.140 

0.778 
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TABLE  I  (Concluded) 
b.  Free-Flight  Data 


M_ 

Re. 

Kn_ 

Tw/T. 

sw 

CD 

CD/CDfnl 

7.65 

4.63 

2.46 

4.69 

2.95 

2.010 

0.828 

7.69 

2.28 

5.02 

4.73 

2.96 

2.340 

0.966 

7.79 

1.11 

10.5 

4.85 

2.96 

2.352 

0.976 

7.59 

2.10 

5.37 

4.62 

2.95 

2.379 

0.979 

7.69 

1.04 

11.00 

4.73 

2.96 

2.447 

1.011 

7.50 

3.82 

2.92 

4.52 

*  2.95 

2.162 

0.878 

7.53 

1.88 

5.99 

4.55 

2.95 

2.300 

0.934 

7.45 

3.40 

3.27 

4.47 

2.95 

2.153 

0.875  ;  . 

7.47 

1.67 

6.66 

4.49 

2.95 
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0.908  1 

7.51 

0.83 

13.40 

4.53 

2.95 
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1.000  ' 

7.51 
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13.20 

4.53 

2.95 
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1.118 

7.37 

3.00 

3.66 

4.38 

2.94 
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0.986 

7.45 

1.45 

7.64 

4.47 

2.95 

•.  2.507 

0.913 

7.29 

2.55 

4.26 
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2.94 
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7.23 

1.30 
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4.23 

2.94 
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2.94 
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4.01 
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6.71 
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2.92 
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6.69 
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15.43 

3.67 
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2.91 

2.440 

1.009 

7.76 

2.22 

5.22 
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8.23 
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7.26 
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2.73 
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5.98 
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2.58 
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4.67 
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